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We investigate in detail the properties of oscillating instanton solutions discussed recently 
in the literature. We find that the solutions with N nodes contain exactly N homogeneous 
negative modes in their spectrum of linear perturbations. The existence of extra negative 
modes for the > 1 solutions suggest that they are not final state physical objects resulting 
from tunneling, but rather unstable intermediate thermal configurations. By contrast, the 
single negative mode for the iV = 1 instanton confirms its interpretation as mediating the 
curved-space tunneling between vacua with equal energy densities. 



I. INTRODUCTION 

Modern developments showed [1, 2] that string theory predicts the existence of a multitude of 
vacua, some of which being stable and some metastable. If this picture turns out to be correct, 
then, via the mechanism of eternal inflation [3], it will have profound consequences for the ultra- 
large-scale structure of the Universe. 

The theory of metastable vacuum decay in flat spacetime was developed long ago [4-9]. It was 
shown that metastable vacuum decay proceeds via the nucleation of bubbles of true vacuum within 
the false vacuum and the subsequent growth of these bubbles. Within the Euclidean approach the 
bubble nucleation process is described by the so-called "bounce" [7], which refers to a classical 
solution of the Euclidean equations of motion with certain boundary conditions. It was shown 
that in the WKB approximation the action of the bounce determines the tunneling rate exponent. 
Furthermore, in flat space-time the pre-exponential factor in the decay rate was calculated by 
taking into account quadratic fluctuations about the classical solution [8, 10-12]. Its value is given 
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by a ratio of the functional determinants of the fluctuation operators corresponding to the bounce 
and the metastable vacuum. It is important to note that there is exactly one negative mode in the 
spectrum of small perturbations about the bounce solution in flat space-time [8, 13]. This single 
negative mode is essential in making the decay picture coherent [14]. 

False vacuum decay with gravity was first investigated by Coleman and De Luccia [15]. It 
was shown that, as in a flat space-time, the Euclidean action of the bounce solution determines 
the leading exponent in the false vacuum decay rate and that the analytic continuation of the 
bounce defines the space-time geometry at the moment of true vacuum bubble materialization. 
The Coleman-De Luccia results were reconsidered and justified by using a technique that explicitly 
accounts for the structure of an initial state of quantum field in the semiclassical calculations of 
the path integral in curved space-time [16] and more recently by a thermal derivation [17]. In 
addition to Coleman-De Luccia bounces, excited multi-bounce solutions have been discussed in 
the literature [18, 19]. It was suggested that in certain regimes these "oscillating" bounces most 
likely will play a significant role in tunneling processes. Instanton solutions mediating tunneling 
between degenerate vacua in curved space were investigated in [19, 21] and excitations thereof were 
studied in [22]. Possible observational tests of cosmological instantons were considered in [23]. 
Gravitational corrections to standard model vacuum decay were studied in [24] and implications of 
recent ATLAS and CMS experimental results about the Higgs boson mass for the stability of the 
electroweak vacuum were investigated in [25]. In spite of much work which has been done since 
Coleman and De Luccia's paper on taking into account gravity in tunneling processes, there still 
remain many open questions, cf. discussions in [26-28]. 

While in fiat space-time finding a negative mode about a bounce is a straightforward task, 
when gravity is taken into account it becomes a more involved problem [29-36]. Using Dirac's 
theory of constrained Hamiltonian systems, it was shown that with the proper reduction procedure 
one finds a single negative mode about Coleman-De Luccia bounces [32, 33]. Furthermore it was 
demonstrated in [35] that the excited multi-bounce solutions possess more then one negative mode 
and consequently they do not contribute directly to the tunneling processes, but rather specify a 
path through configuration space that connects thermally excited horizon volume configurations 
(for a discussion of this point see [17]). 

In contrast to bounces, instanton solutions in flat spacetime typically only have zero modes, and 
no negative modes. As such, they describe the quantum mixing between equal energy states of the 
system (see e.g. the book [37]), rather than describing the decay of one vacuum into another vacuum 
of equal energy density. The aim of the present investigation is to better understand what happens 
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when gravity is included. For this purpose, we will study in detail the properties of the oscillating 
instanton solutions discussed in the literature recently [21, 22]. As we will show, in the presence 
of gravity the standard instantons, which interpolate monotonically between one vacuum and the 
other, possess exactly one negative mode, and they do indeed describe a tunneling process between 
degenerate vacua. By contrast, oscillating instantons contain additional negative modes which 
invalidate their interpretation as contributing to a decay process. Along the way, we will clarify a 
number of issues that arise in the limit where the minima of the potential become degenerate in 
height, such as the proper definition of the thin-wall approximation and the interpretation of the 
existence of the bubble wall itself. 

The rest of paper is organized as follows: in the next Sect, we recall the main properties of 
instantons and bounces using simple quantum mechanical examples. In Sect. Ill we review details 
of the process of tunneling with gravity and we also review gravitational instanton solutions and 
their properties. In Sect. IV we present our analysis of linear perturbations and their spectrum 
about instanton solutions. Sect. V contains concluding remarks. 

II. EUCLIDEAN SOLUTIONS IN FLAT SPACE-TIME: INSTANTONS VS BOUNCES 

Let us consider the basic example of one dimensional field theory (quantum mechanics) and 
assume that the potential V{x) has a double well shape with the degenerate vacua at x = ±a, see 
Figure 1. Classically, the ground state of a particle at rest in one of the minima is a stable state. 
Quantum mechanically it is well known that the existence of instantons, i.e. finite action Euclidean 
solutions interpolating between —a and a , leads to a splitting of the ground state energy (see e.g. 
the detailed discussion in [9], which we briefiy follow in the present Sect.). 

In the semiclassical approximation, summing multi-instanton configurations one finds for the 
two lowest-lying energy eigenstates 



where the exponent S = S[x] is the Euclidean action of the instanton solution and the pre- 
exponential factor 



is calculated by taking the Gaussian integral over the quadratic action describing linear pertur- 
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Figure 1: The left panel exhibits the shape of a double well potential with two degenerate minima located 
at a; = ±a. The blue curve represents the instanton solution, whose profile is depicted in the right panel. 

bations about the instanton. There is one translational zero mode in the spectrum which needs 
special care and det' means that this zero mode is treated separately. Integration over it gives the 
normalization factor in (2). Thus, tunneling effects remove the degeneracy and shift the ground 
state energy. However, the correction to the energy levels is real, and two classically stable states, 
localized around either minimum, are promoted to two non-localized, stationary quantum states. 

When the minima of the potential have different energy, the semi-classical analysis yields quali- 
tatively different results. Now there is one vacuum with higher energy, which is a false (metastable) 
vacuum, and another with lower energy which constitutes the true (stable) vacuum. In this case, 
the bounce solution describes the decay of the false vacuum which proceeds via the nucleation 
of true vacuum bubbles in the false vacuum. In the semi-classical approximation, summing the 
contributions of multi-bounce configurations one finds the following correction to the energy of the 
false vacuum 



[i+om 



(3) 



where 





It is remarkable that there exists exactly one negative mode in the spectrum of linear perturbations 
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about the bounce solution. This can be inferred from the fact that the (translational) zero energy 
wave function -00 ~ of the corresponding Schrodinger equation has a node. The negative mode 
impHes that the correction (4) to the false vacuum energy is purely imaginary and thus we actually 
have a decay process. So, the false vacuum decay probability per unit time is given by 

r = -2imi;/.= (A) 

Analyzing metastable vacuum decay processes in the 1988 NPB article "Quantum Tunneling and 
Negative Eigenvalues" [14] Coleman arrives at the strong conclusion: "There may exist solutions 
in other ways like bounces and which have more than one negative eigenvalue, but, even if they do 
exist, they have nothing to do with tunneling." So, determining the number of negative modes in 
the spectrum of perturbations is of great importance in finding the proper physical interpretation 
of any given solution. 

III. TUNNELING WITH GRAVITY 

Let us now consider a self-interacting scalar field theory minimally coupled to Einstein gravity 
in four dimensions: 

S = I d^x^g (^^R - ^V^^V^^ - V{ip)^ , (6) 

where n = SnG is the reduced Newton's constant. We consider potentials V{(p) which are bounded 
from below and admit two (possibly degenerate) vacua at ip = (p±, separated by a barrier whose 
top is located at ip = iptop (see Figure 2). We denote V± = V{lp±) and conventionally choose 
V- <V+- In other words, when the vacua are non-degenerate, (p- labels the position of the true 
vacuum in field space. 

A. Field Equations 

In the semi-classical limit, the rate for the tunneling process is again described by solutions of 
the Euclidean field equations with proper boundary conditions. In flat spacetime it was shown [38] 
that the 0(4)-invariant solution has least Euclidean action, and provides the dominant contribution 
to the decay process. An analogous proof does not exist in the case where gravity is included - 
rather the 0(4) invariance is still a conjecture [15] that we will assume to be true. Thus, the 
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Figure 2: The left panel shows the shape of a double well potential with two non-degenerate vacua. The 
gray curve represents the inverse potential, driving the solution of the Euclidean field equations. In the 
right panel, the flat space bounce appears as the separating solution between undershooting (red curve) and 
overshooting (green curve). 



Euclidean metric and scalar field Ansatze can be written as 



(7) 
(8) 



The field equations read 



J'2 



+ V) , 



1 + 
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(9) 
(10) 

(11) 



where a prime denotes differentiation w.r.t. Euclidean proper time t] and V^^ = Taking the 
origin p = to be conventionally located at r/ = 0, regularity imposes the following conditions: 



p'ir] = 0) = 1 
if'iv = 0) = 



(12) 
(13) 



Together with (10), these equations can be thought of as describing the one-dimensional motion 
of a unit mass particle which starts at rest at rj = and is subject to an inverted potential 
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—V{ip). According to the sign of p' /p, a friction/anti-friction term is also present, which depends 
non-linear ly on the motion of the particle via (11). 

When ipo = ip{r} = 0) is a stationary point of V , the corresponding solutions are maximally 
symmetric with a constant scalar field profile = ipo^. Depending on the sign of Vq = V[lpq), one 
obtains 

Vb > : Euclidean dS {S^) p = H'^ sin [Hi]), = ^ , (14) 
Vq = : Euclidean space (£'4) p = r] , (15) 



Vb < : Euclidean AdS p = ^ sinh ^^r/j . (16) 

Solutions with a non-trivial profile for the scalar field and which interpolate between the true 
vacuum and the false vacuum, i.e. for which we impose the boundary conditions 

^0 = ^{r} = 0) ~ , 
Lp=ip{r] = fi) ~ if^ , 

also exist. If they contain a single negative eigenmode in their fluctuation spectrum, they describe 
the decay of one vacuum into the other. The value of f/ can be either finite, as for the maximally 
symmetric solution: 

r? = ttH-' , (17) 

or infinite, as in the other cases Vq < 0. We refer to these two types of solutions as compact and 
non-compact bounces respectively. 

B. Bounces Connecting dS Vacua 

Let us focus on the case where both vacua are of the de Sitter type, and the potential is bounded 
by a positive constant, as in Figure 2. In [19] it was shown that the corresponding theories can 
admit, depending on the details of the potential, several compact bounce solutions. These can 
be labeled by the number of scalar field oscillations between the two vacua. Relying on analytic 
arguments, the number of such solutions was determined to be equal to the largest integer N such 



^ By taking derivatives of (10) one can easily see that all the derivatives of tp vanish at = 0. 
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that (see also [20]) 

iV(iV + 3)<%ip^, Hl^^^^^^. (18) 

-"top 

We may obtain a qualitative understanding of the existence of such solutions by analyzing the field 
equations. First, one can easily prove that, V being bounded from below by a positive constant, all 
solutions to the field equations that are regular at = (see (12) and (13)) are compact. Indeed: 

p" = - Y {^'' + V)< -Hi,^ p, Hi,^ = ^min(y) > , (19) 
p(0) = . (20) 

One can view (19) as the equation of motion of an oscillator with a time-dependent frequency, 
bounded from below by a positive constant. It is easy to show that, given the initial conditions 
(20), p vanishes before the timescale set by Hmin- 

p{f)) = : f)< ttH-I. (21) 

However, for a generic value of ipo, one obtains a singular solution. Indeed, when ij approaches fy, 
the coefficient p' /p in the scalar field equation becomes large and negative. For this reason, except 
for a discrete set of values for ipo, the scalar field diverges (positive or negative) when rj approaches 
f/. If LfQ equals one of those special values, the scalar field approaches a constant (pQ when r] ^ r). 

Now, let us try to understand how these particular values of ipQ emerge from the field equations. 
Like in the flat space case, setting (pQ arbitrarily close to the true vacuum one can make sure 
that the corresponding solution for ip overshoots monotonically. However, this does not require the 
two vacua to be non-degenerate anymore. Indeed, if (pQ is sufficiently close to the coefficient 
p' / p will start taking large negative values at r/ ~ ttHZ^, when the scalar field still sits near the 
true vacuum. The combined action of the potential and the anti-friction term will then kick the 
scalar field out of the true vacuum and make it reach infinity at rj = f]. When ipQ — is increased, 
the scalar field "particle" may lose, during the phase in which p' /p is positive, enough energy 
to be trapped and oscillate for a while around the minimum of —V located at (/? = ftop, before 
being kicked out again by the anti-friction at r/ ~ H'^. The period of such oscillations is roughly 
equal to the inverse curvature of the potential at (ptop- The maximum number of oscillations is 
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approximately set by the ratio between these two periods: 



(when iV > 1) . 



(22) 



H 



Regular oscillating bounce solutions appear as the separating elements between classes of solutions 
which have a different number of oscillations before the overshooting: this rough estimate is a 
qualitative version of the rigorous result (18). 



In the limiting case V+ = a Euclidean solution in flat spacetime does not exist. However, 
when gravity is included, instanton solutions between degenerate vacua become possible. Heuris- 
tically, this can be understood as follows: bounce solutions between dS vacua never start and end 
precisely at a minimum of the potential, but rather they start and end a little higher up on the 
potential. This can be interpreted by taking into account the fact that dS space has a non-zero 
temperature. Then, in a path integral formalism the preferred path connecting the original false 
vacuum to a bubble of true vacuum immersed in a region of false vacuum involves first a thermal 
excitation part of the way up the potential, followed by the bounce solution bringing the field over 
the hill. Because of this, the form of the bounce solution itself is insensitive to the behavior of the 
potential outside of the field range that it interpolates along. In other words, the precise height 
of the potential minimum is not directly relevant in determining the shape of the bounce solution, 
and tunneling ought to be able to take place irrespective of whether the minima are at slightly 
different or equal heights. 

As an example of gravitational instanton solutions, we will consider the symmetric potential 
studied in the work of Lee et al. [22], where such instanton solutions were presented: 



The two degenerate vacua are located at ip± = ±1, where V± = Vq. To reproduce the results 
presented in [22], we take Vq = 0.5, k = 0.04. From (18) we find = 6, suggesting that solutions 
with up to six interpolations of the scalar field should exist. The results of the numerical solutions 
of the field equations are represented in Figure 3. The solutions that separate classes of profiles 
with the same number N of oscillations are the regular bounce solutions already described in [22]. 



C. Instantons in a Double Well Potential with Degenerate Minima 
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Figure 3: In the top panel, number of oscillations as a function of the starting value of the scalar field 
fo {Vq — 0.5, K — 0.04): the values separating different classes of solutions correspond to regular instanton 
solutions. In the center panels, profile of the scalar field and p for under/overshooting around the = 1 
instanton. In the bottom panels, same for N — 2. As (po increases, the red profile of the center panels 
evolves towards the green, then blue and red profiles of the bottom panels. 
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In particular, the object separating singular solutions with N oscillations and + 1 oscillations 
is a regular bounce with oscillations. Such regular solutions appear to be always symmetric 
(A^ even) or antisymmetric (A^ odd) under rj ^ fj — rj. This result looks natural, considering the 
symmetry of the potential, but is not predicted by any mathematical proof that we are aware of. 

We would like to stress that the existence of these solutions is not immediately apparent from 
the original formalism developed by Coleman and De Luccia. Much of their work assumes that the 
energy difference between the two minima is small, yet non-zero, in such a way that the bounce 
solution can be effectively thought of as a bubble of true vacuum separated by a thin wall from the 
false vacuum. In the case of compact solutions, a new scale enters the picture, namely the Hubble 
scale, and one must be careful in defining the thin wall limit. 

In flat space the thin wall limit is defined by the condition 



where is a constant, of dimensions (length)"^, which characterizes the potential barrier. When 
the difference in energy density between the two vacua goes to zero, the size of the bubble diverges 
while the thickness of its wall remains approximately constant. In de Sitter space the situation 
is radically different: as we recalled above, bounce solutions survive when the two vacua become 
degenerate. If we require the wall to be thin compared to the bubble size we need, in particular, 
its thickness to be small compared to the size of the whole bounce geometry, which is set by 
Hq^. With this assumption, it is easy to deduce from the overshooting-undershooting argument 
that the wall should be located close to the maximum of p{ri), where the coefficient p' / p vanishes. 
Therefore, in the case of degenerate vacua, the thin- wall regime is attained when the following 
inequality holds: 



In other words, the radius of the bubble approaches the Hubble radius as the energy difference 
between the two vacua tends to zero (see Figure 4). 

Also, in the standard treatment of Coleman and De Luccia, the existence of a bubble wall is 
physically justified by the fact that it carries the difference in energy between the true and false 
vacua. In a sense, the overall energy is conserved in the tunneling process, and the energy that is 
lost by forming a bubble of true vacuum is compensated by the energy of the wall. In our case, the 
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Figure 4: The left panel shows the profiles of the first instanton for the potential in (23) with Vq = 0.5, 
for the values k = 0.4 (blue line), k — 0.09 (purple line) and k — 0.02 (beige line). In the right panel, the 
corresponding profiles of the p function are plotted. As k — 0, the wall position scales as the Hubble radius 
Hq^ cx while its thickness remains approximately constant. 

interior and the exterior of the bubble consist of vacua which have the same energy density, and 
so one may wonder how a wall, carrying tensile energy, is possible at all. The resolution is that, 
in the limit that the vacua become degenerate, a new effect becomes apparent^. In fact, what we 
find is that the presence of the wall is compensated by the removal of a section of the Euclidean 
four sphere. Let us elaborate on this statement: the wall is approximately located at the equator 
of the four sphere, where 

P ^ H,^ ■ (26) 

Let r]A < rjB such that the wall is completely included in the region rj^ < ij < rjB and effectively 
y ~ Vq outside of this interval. We then have 

p ~ Hq^ sin {Hot]) : i] < 7]a , (27) 
p ~ Hq^ sin [Hoir] + Ai])) : t] > i]b ■ (28) 

As we argued above, the thin wall should be located at the equator of the four sphere, where 



Presumably this effect is also present, though likely increasingly subdominant, when the vacua are taken to be at 
heights that differ more and more. 
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Figure 5: Left panel: the scalar field profile for the first, iV = 1, instanton solution with Vq — 1 and 
K = 0.0125. Right panel: the red, solid line represents the corresponding profile of p{ri). The black, dashed 
line represents the function po describing the vacuum de Sitter solution ip — ±1. The black, dotted line 
corresponds to the shifted de Sitter solution (28) with Ar] estimated according to (34). 

p ~ Hq^. If we take rjB — t]a much smaller than H^^ ^ from the field equations we get 



(29) 



We can now estimate the difference 



A^b) - p'oiV- 



•B 



(30) 



As in the standard thin wall limit, the friction/anti-friction term in the scalar field equation can 
usually be neglected at the wall: 



1 



- y ~ -Vq . 



(31) 



Comparing with (28) we find that 



ij— 1 rriB 
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(32) 
(33) 
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that is 



(34) 



This formula gives the width of the shce of four-sphere that is removed due to the presence of 
the wah. Numerical results for the potential (23) indicate a good agreement with this analytic 
consideration (see Figure 5 for an illustration). Thus, loosely speaking, the wall is created by the 
conversion of a slice of spacetime into scalar field gradient energy. 



As discussed above, the study of the spectrum of perturbations about a given Euclidean solution 
is important in order to determine whether the solution actually contributes to the decay rate of 
a given vacuum. Indeed, the presence of a single perturbation mode with a negative eigenvalue is 
required for this interpretation [14]. The negative mode of the Coleman - De Luccia bounce was first 
found in [32, 33] and the fluctuation spectrum about oscillating bounces was studied in [35, 36]. 
In the case of bounce solutions (connecting two non- degenerate dS vacua), numerical evidence 
showed that the A^**^ oscillating bounce is characterized by N homogeneous, 0(4)-symmetric, 
negative modes. For A > 1, there also exist additional inhomogeneous negative modes with non- 
zero angular quantum number £ > 0. 

In their study of oscillating instanton solutions connecting degenerate vacua, the authors of 
[21, 22], probably based on the properties of instantons in quantum mechanics, argued that such 
solutions should not admit negative modes. In what follows we address this question numerically 
and show that this expectation is not fulfilled. 

The quadratic action for linear 0(4)-symmetric perturbations about Euclidean solutions in a 
self-interacting scalar field theory coupled to gravity was derived in [32, 35], and reads: 



IV. SPECTRUM OF LINEAR PERTURBATIONS 




(35) 



with the potential 
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where / = ^'-Q-6ip represents the fluctuation of a scalar field and Q = 1 ^ . The Schrodinger 

equation diagonalizing the above quadratic actions reads 

-f" + U[p{v),^irj)]f = Ef. (37) 

Now our aim is to study the number of negative modes in the spectrum of (37) for different oscil- 
lating instantons. A simple method of counting states with negative eigenvalues of the Schrodinger 
equation consists in investigating the zero energy wavefunction of this equation. According to well 
known theorems [39] , the number of nodes of the zero energy wavefunction then counts the number 
of bound states in a given potential (heuristically, this is because the number of nodes points to the 
existence of a corresponding number of lower energy wavefunctions with successively fewer nodes). 

Solving numerically the background equations for the oscillating bounce solution in the potential 
(23) and the perturbation equations with appropriate initial conditions (see [35] for details) for wide 
class of investigated solutions we find that the instanton with nodes has exactly A^ homogeneous 
(0(4)-symmetric) negative modes. 

First, we considered the bounce solutions connecting two degenerate dS vacua. Our findings 
are illustrated in Figure 6. For definiteness, we have taken the same parameter values as those 
considered in [22], namely Vq = 0.5 and k = 0.04, for which six oscillating bounce solutions are 
possible. The profile of the potential U (??) for the first three solutions is represented in Figure 6. 
It is interesting to note that around the A^-oscillating solution, the potential U admits A^ negative 
minima. The corresponding zero energy wavefunction shows exactly A^ nodes^ (right panels of 
Figure 6). 

Our findings suggest that the results obtained in [35, 36] do not undergo significant changes when 
the energy difference between the two vacua vanishes. This can be made explicit by numerically 
studying the degenerate limit. To do this, we add a linear term to the potential (23): 

V{^) = ^{^^-if + Vo + e^, (38) 
V+-V- = 2e + 0{e^) , (39) 

and investigate the limit e — t- 0. The results of this analysis for a particular value of Vq and k are 

^ Generally speaking, the precise shape of modes and the corresponding eigenvalues depend on the choice of the 
weight function used to define the eigenvalue problem. However, the number of negative energy eigenstates is 
independent of this choice. Therefore, adopting e.g. the weight function of [36], one obtains again the same 
number of negative modes, though in general they differ from the ones computed according to (37). 
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Figure 6: In the left panels, profile of the scalar field for the first three oscillating instantons with N = 1, 2, 3 
and Vq = 0.5 and k = 0.04. In the central panels, potential for 0(4)-synimetric perturbations. In the right 
panels, zero mode wavefunction (dotted line) and negative modes (solid lines): the normalization of the 
wavefunctions is not imposed, so the overall scale of the vertical axis is irrelevant. For iV = 4, 5, 6 we found 
analogous results. 

summarized in Figure 7. We chose the parameters in such a way that the thin wall approximation 
applies to the bounce/instanton solutions. When e — J- 0, the wall moves towards the equator of 
the 4-sphere. However, in the same limit, the negative eigenvalue remains finite, and reaches the 
instanton value. Finally, the wall position can be predicted by the Coleman-de Luccia treatment, 
keeping in mind that the vacuum energy is now finite. 
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Figure 7: Left panel: position of the wall of the = 1 bounce normalized to the size fj of the de Sitter 
sphere, as a function of the energy difference between the two vacua, for Vq — 0.5 and k = 0.4 (dots). t^^oO 
was conventionally defined requiring ip{r]waii) = ftop- The solid line is the wall position as predicted by the 
Coleman~de Luccia treatment. When the two vacua are degenerate, the wall is located at the equator of 
the sphere {riu,aii/v = 1/2)- Right panel: eigenvalue of the 0(4) symmetric negative mode associated to the 
= 1 bounce. 

We also extended our study to bounce solutions connecting Minkowski and AdS degenerate 
vacua, which were also discussed in [22]. As it stands, our formalism only appUes to the cases 
where the function Q is positive. The cases that we investigated are similar in spirit to the 
instanton solutions connecting dS vacua, in that the instantons themselves are compact due to the 
positive value of the potential at ^ptop- It is then not surprising that similar results are obtained for 
these cases, and again A^-oscillating solutions possess N negative modes. It would be interesting to 
see if one can extend this class of solutions to the case where they start and end at zero or negative 
values of the potential. We leave this question for future work. 



V. CONCLUSIONS 

In the present paper we were interested in the properties of oscillating instanton solutions 
in a scalar-gravity theory. We investigated hnear perturbations about instantons and studied 
numerically the corresponding Schrodinger equation. Our results imply that instanton solutions 
with N nodes admit exactly N homogeneous negative modes. The existence of additional negative 
modes for the oscillating instantons with more than one node gives us ground to discard them as 
physical final configurations. They should rather be interpreted as unstable intermediate thermal 
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configurations interpolating between the basic instanton and Hawking-Moss solutions, similarly to 
the oscillating bounces [17]. In contrast to this, the existence of a single negative mode for the 
= 1 instantons supports their decay interpretation. 

We would like to stress that the tunneling process between vacua of equal or similar energy 
densities may not be of solely academic interest: recent studies of eternal inflation have in fact made 
implicit use of such tunneling processes and have even shown that, under certain circumstances, 
these similar-height tunneling processes can be the dominant tunneling processes in the multiverse 
[40, 41]. In an equal-height decay process, the vacuum energy inside of the nucleated bubbles 
will be identical to that in the starting vacuum, yet other physical quantities will of course be 
different in general. In a string theoretic context, the value of the scalar field (p can determine 
certain properties, such as coupling constants, of the low energy physics, and thus there will in 
general be a clear physical distinction between the two vacua. Also, if one considers a more realistic 
model where other fields are coupled to the tunneling scalar field one can observe other interesting 
phenomena like particle creation during tunneling and even more dramatic changes, e.g. if fermions 
have Yukawa couplings to the tunneling field [42, 43]. It is interesting to observe that the rates 
of tunneling in between the two vacua will be the same, regardless of whether one tunnels from 
the first to the second, or from the second to the first vacuum. This is because the tunneling is 
mediated via the same instanton solution, and the background subtraction that one performs in 
calculating the tunneling rate is identical in both cases, as it depends only on the vacuum energy. 
Thus, such a system will result in a fractal spacetime structure with the overall volume being 
divided equally between the two vacua in the late future limit. 

We conclude with a few notes on future directions: the oscillating instantons which we investi- 
gated in present paper exist only due to gravity. There are no such solutions in (3+l)-dimensional 
theory in fiat spacetime. On the other hand in scalar field theory there are so-called Fubini instan- 
tons [44] describing tunneling without barrier. Recently gravitating versions of Fubini instantons 
were investigated in [45]. Since in this case the solutions exist in fiat as well as in curved spacetime, 
it will be interesting to apply a similar analysis to the Fubini instantons and see what the effects 
of gravity are in that case. 

On a more technical side, we point out that in the present paper we studied a wide class of in- 
stantons with the function Q entering the effective fiuctuation potential being positive everywhere. 
In the situation where Q becomes negative somewhere along the trajectory, the perturbation po- 
tential becomes singular and the validity of its derivation needs a more careful analysis. We leave 
these subtle questions for further investigation. 
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